Experimental demonstration of measurement-based quantum computation in 

correlation space 



Wei-Bo Gao*,^ Xing-Can Yao*/ Jian-Ming Cai,^'^ He Lu,^ Ping Xu,^ 

Tao Yang/ Yu-Ao Chen,^ Zeng-Bing Chen,^ and Jian-Wei Pan^'"* 

'Hefei National Laboratory for Physical Sciences at Microscale and Department of Modem Physics, 
University of Science and Technology of China, Hefei, Anhui 230026, China 
'^Institut fiir Quantenoptik und Quantenmformation, 
Osterreichische Akademie der Wissenschaften, Technikerstrafie 21A, A-6020 Innsbruck, Austria 
'^Institut fiir theoretische Physik, Universitdt Innsbruck, Technikerstrafie 25, A-6020 Innsbruck, Austria 
Physikalisches Institut, Ruprecht-Karls-Universitdt Heidelberg, Philosophenweg 12, 69120 Heidelberg, Germany 

(Dated: April 28, 2010) 

The paradigm of measurement-based quantum computation opens new experimental avenues to 
realize a quantum computer and deepens the understanding of quantum physics [JQ • For years, the 
entanglement properties of cluster states have been considered critical for universal measurement- 
based quantum computation Surprisingly, a novel framework namely quantum computation 
in correlation space [lOHl^ has opened new routes to construct quantum states possessing entangle- 
ment properties different from cluster states while still retaining the universality for measurement- 
based quantum computing. The scheme not only offers more flexibility to prepare universal resources 
for quantum computation, but also provides a different perspective to study the fundamental prob- 
lem regarding what are the essential features responsible for the speedup of quantum computers 
over classical devices [l3 - [l6l [. Here we report an experimental realization of every building block of 
the model of measurement-based quantum computation in correlation space. In the experiment, we 
prepare a four-qubit and a six-qubit state, which are proved different from cluster states through 
two-point correlation functions and the single site entropy of the qubits. With such resources, we 
have demonstrated a universal set of single-qubit rotations, two-qubit entangling gates and fur- 
ther Deutsch's algorithm. Besides being of fundamental interest, our experiment proves in-principle 
the feasibility of universal measurement-based quantum computation without cluster states, which 
represents a new approach towards the realization of a quantum computer. 



Measurement-based quantum computation (MQC) 
with cluster states, also called one-way quantum compu- 
tation, has generated enormous interest in the quantum 
information community since its discovery in 2001 [ll-Q. 
In this computational model, two steps are required (i) 
preparing an algorithm-independent universal resource 
state and (ii) performing single-qubit measurements with 
classical feed-forward of their outcomes. Cluster states 
Q have been the only known universal resources for years 
and so far all the reported experimental demonstration of 
one-way quantum computation has been based on cluster 
states U^al • Surprisingly, it was found recently by Gross 
and Eisert that, in the framework of MQC in correlation 
space many singular entanglement properties of cluster 
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states can be relaxed for a universal resource 



The discovery has greatly enriched the resources for 
the implementation of MQC and provided a different 
perspective to study the fundamental problem, that is, 
what are the essential features of quantum states respon- 
sible for the speedup of quantum computers over classical 
devices |13l - [la |? In MQC, the computational power is 
embedded in the universal resource states, entanglement 
properties of which can thus provide insights into the es- 
sential role of entanglement in the speedup of quantum 
computation. Entanglement properties of cluster states 
are usually regarded as being paramount and favorable 
for quantum computation. For example, in cluster states. 



two generic qubits are uncorrelated which makes it 
possible to logically break down a large scale MQC into 
small components, and every particle is maximally en- 
tangled with the rest of the state in order to guaran- 
tee deterministic operations. Nevertheless, in contrast to 
cluster states, universal resource states with arbitrarily 
small local entanglement do exist, and a no n- vanishing 
correlation length as a matter of fact is not an obsta- 
cle for universal MQC [TMl- These striking facts offer 
more flexibility to tailor universal resource states for var- 
ious kinds of physical systerns, e.g. cold atoms and polar 
molecules in optical lattices [17| , or even many-body sys- 
tems in condensed matter physics jl8|. 



Here we report an experimental realization of MQC in 
correlation space. First, we prepare a 4-photon 4-qubit 
state entangled with polarized photons and a 4-photon 
6-qubit state entangled in both the polarization and spa- 
tial modes. We show that they have entanglement prop- 
erties different from cluster states by measuring the two- 
point correlation functions and single-site entropy in the 
states, although they still satisfy the same entanglement 
criteria for universality [3, [3| • Based on the generated 
four-qubit states, single qubit rotations together with the 
working principle of compensating the measurement ran- 
domness are demonstrated. Furthermore, a two-qubit 
entangling gate is implemented based on the six-qubit 
state. Our experiment shows that quantum states differ- 



2 



ent from cluster states can also serve as promising candi- 
dates for one-way quantum computation, and fundamen- 
tally not all entanglement properties of cluster states are 
indispensable for quantum computation. 

General universal one dimensional (ID) computational 
wires are expressed as matrix product states (MPS) [l9l - 
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where ^[|si)] are single-qubit operators, \l) and \r) rep- 
resent left and right boundary vectors in the so-called 
correlation space, which is the auxiliary two-dimensional 
vector space the matrices ^[|si)] act on. A local projec- 
tive measurement on site i with the outcome will in- 
duce the action of the operator = X)(</'iki)^[|si)] 

Si 

in the correlation space, which can serve as the initializa- 
tion, single qubit rotation and readout of a logical qubit 
Quantum ID wires can be coupled to form a 2D 
resource state for the implementation of entangling gates 

il- 

In our experiment, the technique of spontaneous para- 
metric down-conversion (SPDC) and linear optical ele- 
ments are used to generate the resource states (see Meth- 
ods for the details). A four-qubit ID MPS (see Fig. la) 
encoded on the polarized photons is used to demonstrate 
single-qubit rotations. Here we take the left and right 
boundary vectors as |/) = |-|-) and |r) = |0), where 
|-|-) = :^(|0) -I- |1)). Moreover, for the first three sites of 
the state, the tensor matrices are [l^ 

A[\H)]=Hcos9 and A[\V)] = HZ sinO. (2) 

Here \H) and \V) denote the horizontal and vertical po- 
larization, which represent the qubits |0) and |1); H 
denotes the Hadamard gate and Z refers to the Pauli 
matrix cr^. For the end site, we adopt a simple alter- 
native by changing the tensor matrix to = H, 
-^[1^)] = Substituting them into equation (1), the 
four-qubit state is explicitly written as 



1^4) = c\H),{c\H) 
+ s\V),{c\H) 



-s\V)Uc\H)\P)+.s\V)\M))u 
s\V)Uc\H)\M)+s\V)\P))u, 

(3) 



where IP) 



-^{\H) + \V)) and \M) 



\v)). 



For simplicity, we denote c ~ cos 9 and s — sin 9 in the 
above equation and the following text. For any 9 that 
cos9 ^ and sin9 ^ 0, such a MPS can serve as a 
quantum wire for single qubit logic gates [l^ . By tuning 
the parameter 9, the entanglement properties can be very 
different. In our experiment we choose the angle 9 as 
7r/6, and thus c = s ~ 1/2. In this state, the 

processing of logical information is implemented in the 
part of qubits with the tensor matrix A\ while the last 
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FIG. 1: a. A four-qubit state for the implementation of single 
qubit rotations. The black dots are associated with the tensor 
matrices A[\H)] = HcosO and A[\V)] = HZsmO, while 
the red dots are associated with B[\H)] = H, B[\V)] = HZ. 
The measurements on the first three qubits will implement 
the quantum circuit on the right, b. A six-qubit state for the 
implementation of a two-qubit entangling gate. The horizon- 
tal lines (1-2-1' and 3-3') represent two computational wires 
corresponding to two logical qubits. To couple them, we pre- 
pare site A as ^{\H) -\- \ V)) and then apply two controUed-Z 
operations on qubits 2-4 and 3-4. The measurements on the 
qubits 1-2-3-4 will implement a C-Phase gate shown on the 
right. The Z ® Z operation will depend on the measurement 
result of qubit 4, and thus they are drawn in dashed squares. 



site serves as readout and has the function of mapping 
the logical information carried by the correlation system 
to the physical qubit as |0)c — > \P)p and |l)c — > |-^)p 

[a. 

To demonstrate the two-qubit entangling gate, a 2D 
resource state is required. As shown in Fig. lb, our 
2D state is constructed by coupling two ID MPS (1-2-1' 
and 3-3') with site 4 via a general scheme to construct 
universal resources from arbitrary computational wires 
(see Appendix). The coupling is implemented by first 
preparing site 4 as | P) , and then applying two controUcd- 
phase operations between qubits 2-4 and 3-4. We use 
the spatial degree of freedom \H') and \V') of photons 
to carry the qubits 1' and 3'. The corresponding ten- 
sor matrices are respectively B[\H')y] = P[|P')3/] = i7, 
B[\V')v] = B[\M')y] = HZ with |P') = -1=(|H') + |y')) 

and \M') = -^{\H') - \V')). The four-photon six-qubit 
state is thus as follows 

l^/ie) = \H)i\y)i2v\v):>,z' + \V)iZ2\tj)i2i'Zs,\y)3,i,' (4) 

where 

1^) = c^\HHH') + cs\HVH') + cs\VHV') ~ s^\VVV'), 
W) = \{c\HH')+s\VV')), (5) 

and Zi represents the Pauli matrix az applied on qubit 
i. 

For the four-qubit state, we have obtained F = 
{4>A \p\ "04) = 0-74 ± 0.01 by estimating the density ma- 
trix of the experimental state. For the six-qubit state, 
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we have also calculated its fidelity F = {iJjq |p| V'e) = 
0.73 ±0.01 by measuring 36 measurement settings in the 
local decomposition of \iIjq) {i/jq] (see Appendix for the 
details). 

Furthermore, we verify that the entanglement proper- 
ties of our resource states are much different from clus- 
ter states by measuring both two-point correlations and 
local entanglement. The two-point correlation between 
qubit i and j is defined as Q^id*)) = (^l^i ® - 
(4'|ai|\E')(\l/|6j|4'), where di and bj arc respectively Pauli 
matrices acting on qubits i and j. In the state Itpi), there 
is non-zero two-point correlation, i.e. Qxx — 0.375, 
while the corresponding one in the cluster state is zero , 
namely Q}^„,, = . max, ^ \ Ql 



max 

aM=X,Y.Z ' 
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we have theoretical values Q^^ 



0. Similarly, in the six- 
0.433, 



xz 



qubit state {ipQ 

Q?zx ~ 0.375. While for the corresponding cluster state, 
the two-point correlations Q^^^^x ^'^'^ Qmax zero 
Experimentally, we find Q].^ = 0.26 ± 0.01 for \ip4); and 
Qli = 0.41 ± 0.02, Q34 ^ 0.32 ± 0.02 for ji/'e)- Although 
the experimental data is not ideal, none of the tested 
two-point correlations equals zero, clearly showing that 
our states arc not cluster states. 

Local entanglement between one single qubit i and the 
other qubits is quantified by the local entropy E{pi) = 
2(1 — Tr£^), where pi is the reduced density matrix of 
qubit i [22|. In a cluster state, every qubit is maximally 
entangled with the rest of the state, and thus the local en- 
tropy always equals 1 0, Q • Instead, in the resource state 
we have theoretically E{pi) ~ E{p3) ~ 0.75 and 
E{p2) = 0.5625. In ji/'e), ^(pi) = ^(^2) = ^(pa) = 0.75 
and E{p4) = 0.9375. Experimentally, in \ip4), the local 
entropy are respectively E{pi) ~ 0.93 ± 0.01; E[p2) = 
0.63 ± 0.01; £;(p3) = 0.88 ± 0.01. In ji/'e), E{pi) = 
0.62 ± 0.02; E{p2) = 0.80 ± 0.02; £:(p3) = 0.67 ± 0.02; 
Eipi) = 0.90±0.02. The experimental data is imperfect; 
however, none of them arc equal to 1, which shows the 
different entanglement properties of our resource states 
from the cluster states. 

To realize an arbitrary SU(2) single-qubit rotation, we 
measure the qubits of the four-qubit state 1-04) in the ba- 
sis B{C) = {|Co), ICi)}, where |Co> ^ s\H 



ictan I \V) 



and ICi) = c\H) — is cot | \ V). For simplicity, we de- 
fine the outcome Vj to be if the measurement re- 
sult is |Co)> and as 1 if the result is We first 
consider the case when all the measurement outcomes 
are 0. This will reduce the success probability of the 
gate for each step of measurement, but it suffices as a 
proof-of-principle to demonstrate the single-qubit gate 
(the strategy to compensate the wrong outcome and im- 
prove the success probability will be demonstrated in 
the following text). In each measurement step, single 
qubit rotation RziC) = exp(— iCo'z/2), followed by a 
Hadamard operation, can be implemented. Based on the 
generated four-qubit state, we perform consecutive mea- 
surements Bi{a), B2{(i), ^3(7) on the physical qubits 



1, 2, 3. By doing so, the input state of the correla- 
tion system |V'm)c = |+) is transformed into |V'oMt)c = 
HRMRx{P)RM)\+), where i?,(/3) = exp(-z/3aj2) 
(See Fig. la). Following the map induced by the ten- 
sor matrix B, the output state of the physical state is 
\i^out)p = Rz{i)Rx{I3)Rz{.q)\+)\o-^h,i->v ■ The experi- 
mental fidelities of six states are shown in Table. I., the 
average of which is 0.86 ±0.01, clearly above the classical 
threshold 2/3 More results of the single-qubit gate 
can be found in the Appendix. 
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TABLE I: The fidelities of tlie output states of the single- 
qubit rotation. The first three qubits of the four-qubit cluster 
state are measured in basis B{a), B{I3), -6(7). \ipout)c and 
\^out)p represent the output states in the correlation space 
and the physical world, respectively. |/o) = (|0) 
and \h) = ^[\Q)^i\l)). 



A distinct feature of MQC in correlation space is the 
strategy to compensate the randomness of measurement 
outcome. In cluster state quantum computation, the 
wrong measurement outcome only induces by-product 
Pauli operators, which can be compensated by using 
feed-forward control of future measurement basis. The 
by-product operators in MQC in correlation space can 
be different from Pauli matrices, and the simple feed- 
forward technique does not work in general any more. 
Fortunately, it has been proved that the introduced er- 
rors can still be efficiently corrected in a bounded number 
of steps with the trial- until-success strategy U 12 1 . 

As a proof-of-principle, we compare the success prob- 
ability of implementing the rotation HRz{a)\+) with a 
two-qubit state IA34) = c \H)^ \P)^+s \V)^ and with 
the four-qubit state \'4>a)- Based on the two-qubit state, 
we measure qubit 1 in the basis B{a). When the out- 
come 7'i = 0, we obtain the desired single-qubit rota- 
tion HRz{a)\+) with a success probability Psia). How- 
ever, when ri = 1, we obtain a rotation HRz{a')\+) 
but with a wrong angle tan(Q;'/2) = (—1/3) * cot(Q;/2). 
Based on the four-qubit state, we first measure qubit 1 
in the basis B{a). When obtaining ri = 0, we mea- 
sure the qubits 2 and 3 in the Z basis, resulting an out- 
put state X''^Z''^HRz{a)\+). When ri = 1, we mea- 
sure the second qubit in the Z basis with the outcome 
r2 and the third qubit in the basis B[{—lY^{a — a')]. 
When r3 = 0, with a probability Ps{o: — a'), we ob- 
tain the output state Z'^^HRz{a)\+). The final result 
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is equivalent to the desired rotation up to a Pauli by- 
product operator. Thus, with two more qubits, we can 
boost the success probabihty of the rotation from Ps{a) 
to Psice) + [1 — Ps(,ci)] • Ps(a — a') (See Appendix for the 
case of more qubits). 
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• Experimental data with compensation 
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FIG. 2: Theoretical and experimental success probability of 
compensating the rotation error with two-qubit states and 
four-qubit states. The blue (red) lines represent the theo- 
retical success probability under the white noise model in the 
two-qubit (four-qubit) case. The blue (red) dots represent the 
experimental values in the two-qubit (four-qubit) case, which 
are achieved by measuring the states for 5s (300s). From the 
figure, we can see that the success probability has been en- 
hanced by using two more qubits. The experimental data is in 
good agreement with the theoretical values when considering 
the effects of white noise. 

The experimental success probability in the four-qubit 
case is achieved by measuring in the appropriate basis 
and then adding the probability of each successful mea- 
surement branch. In the analysis of the experimental 
data, wc calculate the corresponding theoretical success 
probability with the states under the white noise model. 
In other words, wc write the two states as: 

h 

P2 =P|A34>(A34| + (1 -P)^^ 

P4=P'|^4>(V'4| + (1-P')^, (6) 

16 

where /2(/4) denotes two-qubit (four-qubit) identity ma- 
trix; p = (4/2 - l)/3; p' = (I6/4 - 1)/15 and /2(/4) is the 
fidelity of the two-qubit (four-qubit) state, that is, 0.90 
(0.73). The theoretical and experimental success proba- 
bilities are shown in Fig. [2l from which we can see that, 
the experimental data is well consistent with the theo- 
retical curve and the success probability has indeed been 
increased with compensation. 

Besides singlc-qubit rotations, a two-qubit entangling 



gate is required to demonstrate universal quantum com- 
puting. Based on the six-qubit state If/'e); we have real- 
ized a two-qubit controlled phase gate. We first measure 
qubit 1 in the basis B{a). When ri = (we only con- 
sider this case in the following text, the other outcome 
just corresponds to a different input), this measurement 
initializes the input state of the correlation system as 
HRz{a)\+)t (81 |+)c, where t and c denote the target and 
control logical qubit corresponding to two computational 
wires (1-2-1' and 3-3'). Then, we measure qubit 2 and 3 
in the basis S(f) = {s\H) + ic\V) ,c\H) - is\V)}. (a) 
Once having the outcome r2 = rs 0, we measure qubit 
4 in the Y basis {\H} + i\V) , \H) - i\V)} and transform 
the logical state carried by the correlation systems from 

\iJ^n) = HRz{a)\+)t(^\+)c to lijout) = imH)iZ(^ZY^- 

CZlV'm), where CZ = I ® \0){0\ + Z ^ is the con- 

trolled phase gate. The tensor matrices corresponding to 
qubits 1' and 3' map the logical output state to the phys- 
ical output state [16[ carried by qubits 1' and 3' as (/ ® 

H) ■ {Zy ® ZyY' ■ (|0)|+) -Ztanf |l)|-))v3'|0-.ffM-.V". 

(b) If r2 7^ or r3 7^ 0, we measure qubit 4 in the Z ba- 
sis and decouple two ID wires. After the measurement, 
a local rotation error will appear on each logical qubit. 
The same as in the single-qubit gate, we can efficiently 
correct them by applying the trial-until-success strategy 
if we have more qubits. 

In the experiment, when r2 ~ 7*3 = 0, wc charac- 
terize the output state of the two-qubit gate by the 
method of state tomography. We collect the experimen- 
tal data for 600s for each of 36 combinations of the mea- 
surement basis {\H') , |T/') , \P') , \M') , \R') , |L')}, where 
\R') = ^i\H')+^\V')) and \L') = ^ QW) - ^ \V')), 
and then estimate the density matrix with the maxi- 
mum likelihood technique. The theoretical and experi- 
mental density matrices of physical output states when 
a = 0, r4 = and a = 7r/3, r4 = are shown in Fig. |31 
The states are in good agreement with the ideal states, 
which can be seen from their fidelities: 0.88 ± 0.02 and 
0.84 ±0.03. 

When r2 7^ or 7^ 0, we have measured the fi- 
delity between the experimental states and the expected 
ones. In the experiment, a is fixed as one of the angles 
of 0, 7r/2, TT, 7r/3, —tt/3. When r2 = 0, rs = 1, the average 
fidehty of the measured output states is 0.88±0.01; when 
r-i = 1,'"3 = 0, the average fidelity is 0.85 ± 0.01; when 
r-z = l,r3 = 1, the average fidelity is 0.87 ±0.01. The re- 
sults are well consistent with the theoretical predictions. 
More data when r2 = r3 = and the detailed fidelities 
when r2 7^ or r3 7^ are shown in the Appendix. 

The realization of multiqubit quantum algorithms is 
the ultimate goal of quantum computation and the exper- 
imental community has made considerable efforts toward 
this aim. Deutsch's algorithm represents an interesting 
instance of demonstrating the power of quantum compu- 
tation over classical computation. It has been experimen- 
tally reported in the quantum circuit model (23 | and the 
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FIG. 3: The density matrices of the output states of the 
entanghng two-qubit gate. The left column shows the real 
and the imaginary part of the experimental output state den- 
sity matrices, while the right column shows the theoretical 
density matrices. In cases a and b, qubit 1 is respectively 
measured in bases B{0) and B{tt/3). All the matrices are 
achieved in the case when the outcome of the measurements 
is ri — r2 = rs = r4 = 0. 



one-way computation model based on cluster states [25| . 
Here we experimentally demonstrate the implementation 
of Deutsch's algorithm based on the six-qubit states IV'e)- 
Deutsch's algorithm, also known as the Deutsch-Jozsa 
algorithm [2^, allows one to distinguish two different 
types of function f{x) implemented by an oracle in a 
black box. The function f{x) with an N-bit binary in- 
put X is constant if it returns the same value (either 
or 1) for all possible inputs and balanced if it returns 
for half of the inputs and 1 for the other half. With 
classical methods, in some cases we have to query this 
oracle as many as 2^~^ + 1 times. However, only one 
query is required in all cases by using quantum computa- 
tion (26j . In the two-qubit version, the applied algorithm 
can be demonstrated as — |x)|y © fi^))-, where 



|x) is the query input qubit and \y) is the ancilla input 
qubit. Preparing as |+)|— ), we can get the output 

states [(-1)^(°)|0) + Finally, we apply 

two Hadamard operations on the two output qubits. If 
/(.t) is constant, the final result will be \H)\V), while 
if it is balanced, the final result will be There- 
fore, now we can determine the types of f{x) from the 
output states. There are two types of constant function, 
namely /i ~ I ® I and /2 ~ I ® ax- Also, there are 
two types of balanced function, namely /a = gcnot and 
/4 = gcnot{I®ax), where gcnot represents a controlled- not 
gate. Since /2 and fi can be achieved from /i and /a by 
single-qubit local operation, in the following we consider 
only the functions /i and /s. 

In the experiment, we use the wire 3-3' of jV'e) to 
carry the query logical qubit and 1-2-1' of l-^g) for the 
ancilla logical qubit. In the following, we consider the 
implementation of Deutsch's algorithm when ri = r2 = 
ra = r4 = 0. The measurement of qubit 1 in the ba- 
sis B{tt) and qubits 2, 3, 4 in the basis B{0) performs 
the transformation Hg (X) Rz{Ti)a- Moreover, considering 
the map of logical qubits from correlation systems into 
the physical qubits the final transformation of the inputs 
will be {Hq ® Ha){Iq ® Rz{TT)a), which implements the 
preparation of ancilla qubit |— ), the constant function 
and the readout of the qubits (See Fig. 4a). Similarly, 
by measuring qubit 1 in the basis B{tt) and qubits 2, 
3, 4 in the basis B{tt/2), we can implement the gate 
{Hq (g) Ha)gcnot{Iq ® Rzi''^)a) On the iuputs, which is the 
case when the function is balanced (See Fig. 4b). Here 
the measurement of qubits 2, 3 and 4, together with the 
detector signal corresponding to ri = ^2 = ra = r4 = 
should be viewed as the entire black box. In the exper- 
iment, the success probability of recognizing the type of 
fix) is as large as 99% ± 1% for /i and 75% ± 2% for 
/a. The non-ideal probability is mainly caused by the 
non-perfect resource state and single-photon interferom- 
eters in the setup. In the Appendix, we have discussed 
in detail the cases when some measurement outcomes are 
not 0. 

In summary, we have realized a proof-of-principle ex- 
perimental demonstration of measurement-based quan- 
tum computation in correlation space. Every building 
block of the scheme and a simple yet interesting al- 
gorithm has been demonstrated, providing evidence of 
promising applications of states different from cluster 
states. These new resource states do not possess many 
apparently necessary entanglement properties of cluster 
states, such as vanishing correlation functions and being 
locally maximally mixed, but nevertheless are universal 
for quantum computation. 

There are many open questions, both theoretical and 
experimental, worth investigating in the future. In 
measurement-based quantum computation with cluster 
states, it is known that if the product of the number 
resolving detector efficiency and the source efficiency is 
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FIG. 4: The theoretical design and the experimental results 
of implementing Deutsch's algorithm. The wire 3-3' repre- 
sents the query logical qubit, and the wire 1-2- 1' represents 
the ancilla logical qubit. In a and b, the measurement pat- 
terns on the physical qubits (left) implement the algorithm 
on the logical qubits (right). The angles a represent the mea- 
surement basis B{a) of each physical qubit. In the circuit, 
the operation (/ d5 H) represents the mapping of the qubits 
from the correlation space to the realistic physical world. The 
final results are carried by the qubits 3' and 1'. The success 
probability of recognizing the constant function is 99% ± 1%, 
while for the balanced function, this probability is 75% ± 2%. 



greater than 2/3, efficient linear optical quantum com- 
putation is possible |27[. With novel resource states for 
quantum computation in correlation space, computation 
will probably require more qubits. However, bonds be- 
tween two particles are easier to create the lower the re- 
quired entanglement [28|. It would be interesting to in- 
vestigate what conditions are required to make the latter 
scheme scalable. For realistic systems, a trade-ofF may 
exist between the effort of preparing a universal resource 
and its efficiency for quantum computation. The feed- 
forward rule also warrants further study, and it would be 
desirable to combine the feed-forward technique with the 
proof-of-principle demonstration. 



METHOD 

Creation of the resource states. 

In our experiment, entangled photons are created by 
using type-II parametric down conversion (29j . Fem- 
tosecond laser pulses (« 200 fs, 76 MHz, 788nm) are 
converted to ultraviolet pulses through a frequency dou- 
bler LiBzO^ (LBO) crystal and the ultraviolet laser pulse 
passes through two nonlinear crystals (BBO), generating 
two pairs of photons in path 1-2 and 3-4. The observed 
two-fold coincident count rate is about 5.4 x 10^/s. 




PBS A, 



1 HH' ■ ' HH' ■ 



^' A, 

FIG. 5: Experimental setup for (a) the generation of the four- 
qubit state 1^)4) and (b) the six-qubit state \4>e)- a. An ultra- 
violet pulse passes through two BBO crystals to create two 
pairs of entangled photons. Then, half-wave plates (HWP) 
and a series of polarization dependent beam splitter cubes 
(PBCs) have been used to create the desired state. Prism Adl 
is used to ensure that the input photons arrive at the PBCs 
at the same time. Furthermore, every output is spectrally 
filtered {AXfwhm ~ 3.2 nm) to ensure good temporal over- 
lap. A combination of HWP, quarter-wave plates (QWP) and 
polarization beam splitter (PBS) has been used to implement 
the measurement setups, b. By exchanging the labels of 1 and 
2 of |?/'4), and letting the photons 1 and 3 enter two PBSs, the 
desired six-qubit state \tp(j) can be prepared. The ultra-stable 
Sagnac-ring interferometers have been used to measure the 
spatial qubits. In the interferometer, the specially designed 
beam-splitter cubes are half PBS-coated and half beam split- 
ter coated and high-precision small-angle prisms are used to 
satisfy fine adjustments of the relatively delay of the two dif- 
ferent paths. 
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First, let's consider the preparation of the four-qubit 
cluster state. By placing half-wave plates, we prepare the 
initial two-photon states as 



and 



l/V2{\H)^\P)2 + \V)i\M)2) 



l/V2{\H):i\P)i + \Vh\M)i). 



(7) 



(8) 



Then, two polarization dependent beam splitter cubes 
[PBCa) with Ty = {s/cY and Th = 1 are placed in paths 
2 and 3, where T„ {Th) represents the transmission prob- 
ability for the vertical (horizontal) polarization. This will 
transform the entangled two-qubit states into 

l/V2[\H)MH) + s\V))2 + \V)MH) s\V))2) (9) 
and 

l/\/2(c|i/)3|P)4 + s|l/)3|M)4). (10) 

Finally, we let photons 1 and 4 enter a series of PBCs. 
According to the method in Ref. [s^ , we can apply a C- 
phase gate between qubit 1 and 4 through a combination 
of PBCs, which includes an overlapping PBCa (Ty = 1/3 
and Th = 1), and a PBCb {T^ = 1 and Th = 1/3) in each 
path of 1 and 4. Note that the function of PBCb in 
path 1 is only to adjust the amplitude of \H)-^ and \V)-y. 
Therefore, by removing the PBCb in path 1 (using only 
the overlapping PBCa, the PBCb in path 4), we can 
prepare exactly the state |V'4)- 

Next, let's consider the preparation of the six-qubit 
state l^/ig)- Based on the created state 1-04), we ex- 
change the labels of the outputs 1 and 2. Then, we let 
each of the photons 1 and 3 enter a PBS (see Fig. 5b). 
Since a PBS transmits H and reflects V polarization, H- 
polarized photons will follow one path and ^-polarized 
photons will follow the other. In this way, two spatial 
qubits are added onto the polarization qubits 1 and 3: 



\HH'),^,, and \V), 



WV) 



where i ~ \,i 



and the levels are denoted as \H') for the first path and 
\V') for the latter path. After these operations, the final 
state is converted to exactly a four-photon six-qubit 
non-cluster state. In the experiment, the spatial qubits 
play the role of reading out the results, and we need 
to measure the spatial qubits in the X, Y bases, which 
requires matching different spatial modes on a common 
BS. We have designed a special crystal combining a PBS 
and a beam splitter, and then we have used Sagnac-ring 
interferometers to construct the required single-photon 
interferometers, which can be stable for almost 10 hours 
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APPENDIX: EXPERIMENTAL DEMONSTRATION OF MEASUREMENT-BASED QUANTUM 

COMPUTATION IN CORRELATION SPACE 

Universal scheme for coupling ID computational wires into a 2D resource 

The scheme we demonstrate in our experiment to couple two ID computational wires and the measurement pattern 
for the implementation of an entangling gate is rather general and applicable to arbitrary computational wires. In 
the following, we reformulate the scheme for general wires, which our experiment provides a specific example of. 

Quantum wires, with tensor matrices in the canonical form Q as A[|0)] = W and = WS{e) where 

S{9) = diag(e~'' , e*^), can be coupled into a 2D resource state which is universal for measurement based quan- 
tum computation. The coupling scheme used in the experiment is depicted in Fig. [1] Qubit 4 is first prepared in 
1+) = 2-^/'^{\Q) + |1)), two controllcd-X gates (in the same basis as the canonical form of the tensor matrices) are 
applied between qubits 2 and 4, qubits 4 and 6. The scheme is universal in the sense that it works for all general 
wires in the canonical form and the applied coupling gates, namely controUed-X gates, are independent of the specific 
quantum wire. To decouple the wires, we only need to measure qubit 4 in the computational basis {|0), If the 
outcome is 0, we have un-done the coupling; otherwise we get the original ID wires, up to the action of Ux on qubits 
2 and 6. 



1 2 3 




• b o 

5 6 7 

FIG. 1: (Color online) Universal coupling scheme for arbitrary quantum wires. Qubit 4 is first prepared in the state 1+), two 
controUed-X gates (in the basis of the canonical form of the tensor matrices) are applied between qubits 2 and 4, qubits 4 and 
6 (red dashed lines). 

To perform an entangling gate, we do the following measurements: (a) Measure qubits 2 and 6 in the basis 
{uo|0) — ui|I), -uilO) +wo|l)}, where uq = 2"-^/^ (cos | — sin |) and ui = 2~^/^(cos | -|-sin |), the outcomes are denoted 
as S2 and sq; (b) If S2 = 1 or sg = 1, measurement of qubit 4 in the computational basis will decouple the two quantum 
wires and the logical qubits, carried by the correlation systems corresponding to the two quantum wires, leaving a 
byproduct operator Ob- We can then restart the procedure after compensating Ob via the trial- until-success strategy 
0, S]; (c) If S2 = sg = 0, we measure qubit 4 in the Y basis {|0) -I- i|I), |0) — i|I)} and implement the controlled 
phase gate with a by-product operator [VF^d ) (g) M^5'(f )] • (cr^ ® a^Y^ ■ {I (g) |0)(0| -f ® |1)(1|) on the correlation 
systems. The trial-until-success strategy to compensate the randomness of measurement is used to guarantee the 
computational efficiency 0, Q when we obtain the wrong measurement outcome in the above protocol. 
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The characterization of the four-qubit state 



To characterize the four-qubit state, we extracted its density matrix by the method of over-complete state to- 
mography. We collected experimental data for 300s for each of the 1296 combinations of the measurement bases 
{\H) ,\V) ,\P) ,\M) ,\R) ,\L)}, where \R) = -^i\H)+i\V)) and \L) = -j^ {\H) - i\V)). With these data, the 
maximum-likelihood technique is used to construct the density matrix of the state. Even though the state may 
look different in Fig. 2, what is essential is the fidelity which characterizes how well the task is performed, 
F = ('04 \p\ -04) = 0.74 ± 0.01. Similar to the method in refs. H El, the error bar of the fidelity is calculated 
by performing a 100 run Monte Carlo simulation of the whole state tomography analysis, with Poissonian noise added 
to each experimental data point in each run. 







FIG. 2: a. The real parts of the measured four-qubit state density matrix, b. The real parts of the expected density matrix, 
c. The imaginary parts of the measured four-qubit density matrix, d. The imaginary parts of the expected density matrix. 



The fidelity of the six-qubit state 



To obtain the fidelity of the six-qubit state |'06)j we use a method similar to using an entanglement witness @, 
We decompose the density matrix |'06)('06| into locally measurable observables, the values of which can be obtained 
by measuring 36 settings. The qubit is carried by the polarization or path degree of freedom of a photon, i.e. 
|0) o \H), \H') and |1) o \V), \V'). The state [V'e) can be given as 

\^P,) = ±i\a) + \b)) (11) 



where 



|a) = 10)4 (c |00) + s |11))33' ® [c 10)2 (c |00) -I- s |ll))n' + s (c |00) - .s |ll))ir]; 

\b) = 11)4 (c|00)-s|ll))33'®[c 10)2 (c|00)+s|ll))n'-s 11)2 (c|00)-.s|ll))n']. (12) 



The decomposition of \ip6) is 



36 

1^6) (06 I =51^'^^. 

i=0 



(13) 
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where 



Ml = /4 ® (c2 100)33, (00| + |11>33' (111) ® (c' |0>2 (0| + |1>2 (1|) ® (c' 100)^ m + Ill>ii' (111) ; (14) 



M2 = Z4® (c^ 100)33, (00| + «^ 111)33' (111) ® ® (c^ |00)ii, (00| + |ll)^^, (il|) ; (15) 

Ms = Zi® QcsXa ® X^^j ® (c^ 10)2 (0| + \\)^ (1|) ® (c^ |00)ii, (00| + (11|) ; (16) 

Mi = h® (c2 100)33, (00| + -"^^ 111)33' (111) ® 11)2 (1| - c2 10)2 (0|) ® icsFi ® Yy; (17) 

M5 = /4®(c2|00)33, (00|+s2|ll)gg, (ll|)®(c2 10)2(01 -s2|i^_^(;LQ^ l^g^^^Jj^^,. (18) 

Afe = -^4 ® QcsFa ® Fg') ® (c' |0)2 (0| + ll)^ (1|) ® (c^ |00)ii, (00| + \n)^^, (11|) ; (19) 

Af7 = /4 ® QcsXa ® X^^j ® CSX2 ® (c^ |00)ii, (00| - |11)^^, (11|) ; (20) 

Ms = h® i^csYs ® ^3'^ ® CSX2 ® (c^ |00)ii, (00| - s2 |ii)^^^ (iiQ . (21) 

Mg = Z4 ® (c^ 100)33, (00| + 111)33' (111) ® ^"^^^^^s ® ^1 ® Vi'; (22) 

Mio - Z4 ® (c^ 100)33, (00| + «^ 111)33' (111) ® ^"^^^^^s ® n ® ^1'; (23) 

Afn ^ Zi® QcsXa ^3.^ ® (c^ |0)2 (0| - |i)_^ (;l|) ^csXi Xy] (24) 

M12 = Z4® QcsFg^yg.^ ®(s2 11)2(11 -0^10)2(01) ®icsXi®Xr; (25) 

Mi3 = Z4 (8 (^csXs ^3'^ ® {s^ 11)2 (1| - 10)2 (0|) ® icsFi (8 Yy- (26) 

Mi4 = Z4 ® QcsFs ® Y^^j ® (c^ 10)2 (0| - 11)2 (1|) ® icsYi (g) y^; (27) 

Mi5 = /4 (g (^"^^^3 ® ^3'^ ® <»Xi® Yy; (28) 

Mi6 = /4 ® QcsFa (8^3'^® ^c2s^^2 ®Xi® Yy, (29) 

Mi7 = /4 (g (^"^^^3 ® ^3') ® ^'^^^^-^2 <E)Yi® Xy; (30) 

Mi8 = /4 (8 (^"=^^3 ®Y3,^® \c^S^Y2 ®Y^® Xy; (31) 

Afi9 = ^4 ® (c2 100)33, (00| - 111)33' (111) ® (c' 10)2 (0| - 11)2 (1|) ® (c' |00)n, (00| + \n),,, (111) ; (32) 

M20 - n ® (c^ 100)33, (00| - 111)33' (111) ® csr2 ® (c^ |00)ii, (00| - |11)^^, (111) ; (33) 

A//21 = X4® (c^ 100)33, (00| - «^ 111)33' (111) ® (c^ 10)2 (0| + 11)2 (1|) ® ^csXi ® Xp; (34) 

Af22 = X4® (c^ 100)33, (00| - 111)33' (111) ® (c^ 10)2 (0| + 11)2 (1|) ® ^csYi ® Yy; (35) 

M23 = IcsYi ®X3® 13' ® (c' 10)2 (0| - 11)2 (1|) ® (c' |00)ii, (00| + \n),y (11|) ; (36) 

M24 = icsr4®l3®^3'®(c'|0)2(0|-s2|i)2(l|)® (c2|00)ii, (00| + s2|ll)^^, (111); (37) 
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Mas = 




111) 


33' (111 


- 100)33, 


(00|) $ 




X2®Xi 


® yi'; 


(38) 


M26 = 




111) 


33' (111 


-c^ 100)33, 


(00|) € 




X2®Yi(. 


3^1'; 


(39) 


M27 = 


^csX4 (2 


)^3 


«) y3' 


?) CSF2 <8) (s 


'|ll>n 


(iii- 




(001) ; 


(40) 


M28 = 




5^3 




3 CSF2 <8) (s 


'|ll>n 


(iii- 




(001) ; 


(41) 


M29 = 




1^3 




)(^' 11)2(1 


-c^ 10)2(01) 






(42) 


A/30 = 




1^3 




) (c2 10)2 (0| 


-.^|1 


)2(1|) 


53 ^csXi 


® Xi,; 


(43) 




1 

-csy4 (g) 


-'^3 


^ y3' « 


) (c2 10)2 (0| 




)2(1|) 


1 


® Xv; 


(44) 


M32 = 


icsF4 ® 






11)2(1 


-c^ 10)2(01) 






(45) 


M33 = 






!8)X3/ $ 


2 










(46) 


M34 = 


^csX^ € 


5i"3 


(SX3, <i 


2 


g) Fi «) 








(47) 


M35 = 


-CSX4 € 


5^3 


® y3, $ 


2 


gXi ® 








(48) 


M36 = 


^csX^ € 


5^3 


® y3, $ 


3 -c^s^X2 
2 




Xy. 






(49) 



To obtain the experimental values of the above operators, we need to measure using the 36 measurement settings 
(4-3-3'-2-l-l'): ZZZZZZ; ZXXZZZ; ZYYZZZ; ZZZZXX; ZXXZXX; ZYYZXX; ZZZZYY; ZXXZYY; 
ZYYZYY] ZZZXZZ- ZXXXZZ- ZYYXZZ- ZZZYXY- ZXXYXY- ZYYYXY; ZZZYYX- ZXXYYX; 
ZYYYYX- XZZZZZ- XZZZXX- XZZZYY- XYXYZZ- XYXXXY- XYZXYX; XXYYZZ; XXYXXY; 
XXYXYX- YYXZZZ; YYXZXX; YYXZYY; YXYZZZ- YXYZXX; YXYZYY- YZZYZZ; YZZXXY; 
YZZXYX. The theoretieal and experimental values of the 36 observables are listed in Fig. [3] 




Ml W2 M3 M4 W5 MB Nf7 MB M9 MID Mil Ml 2 M13 M14 M15 Ml 6 Ml T M1B MlS M20 ^(21 M32 M23 WSfl M35 M26 M2? M2a M29 M3l M32 MSS M3d M35 MSB 

Observable Observable 



FIG. 3: The theoretical and experimental values of the required 36 observables. a, The required 18 observables to obtain the 
value of \a){a\ + \b){b\\ b, the required 18 observables to obtain the value of \a){h\ + |b)(a|, where \a) and \b) are states in Eq. 
(2). 
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Additional experimental results of the single-qubit gate 



Based on the generated four-qubit state, an arbitrary SU(2) single-qubit rotation can be implemented by the 
consecutive measurements Bi{a), B2{f3), -63(7) on the physical qubits 1, 2, 3, transforming the correlation system 
from the state \ipin)c = |+) to \'ipout)c = HRz{j)Rx{P)Rz{a)\+) , with the corresponding physical output state being 
\'4'out)p = Rz{'y)Rx{P)Rz{ct)\+)\o^H,i^v ■ In Table I, we give additional experimental results of the fidelities of the 
output states. 



a P -y \i'out)c \tpont)p fidelity 

TT 7r/3 ^10) + ^|1) i|P) + ^|M> 0.78 ±0.01 

-7r/3 ^|0) + ||1) ^ IP) + II Af) 0.88 ±0.01 

7r/2 27r/3 7r/2 ^|0)±1|1) ^ |P) ± i |M) 0.91 ± 0.02 

tt/2 tt/3 tv/2 i|0)±^|l) l\P) + ^\M) 0.84 ± 0.02 

7r/2 27r/3 - ^|+) i|-l/> - 0.70 ± 0.02 

n/2 n/3 ^h) - ^1+) ^IV) - h^) 0.74 ± 0.02 

TABLE II: The fidelities of the output states of the single-qubit rotation. The first three qubits of the four-qubit cluster state 
are respectively measured in bases B{a), B{I3), -6(7). \ipout)c and \ipout)p represent the output states in the correlation space 
and the physical world, respectively. In each case, the fidelity is achieved by measuring in 20 minutes. 



Trial-until-success strategy to compensate the randomness of measurement outcome 

The implementation of an arbitrary SU(2) single-qubit rotation can be achieved by measuring the qubits, with 
the tensor matrices = Hcos9 and ~ HZsind, in the following basis B{a) — {|ao),|ai)}, where 

|ao) = sin9\H) +i cos ^? tan ^ |^) and |q;i) = cos 6* \ H) ~ i sin cot f If the outcome r = 0, we realize the desired 
operation HRz{a) with the success probability as 

, , sin^26' sin^26' , , 

= 2(1 -cos 20 cos a) - 2(l±|cos20|) = ^^^^ 

Thus, once sin(20) ^ 0, the success probability Ps{c() is always lower bounded by a positive constant (i.e. independent 
on the angle a). If we get the wrong outcome r = 1, with the probability 1 — ps{a), we actually implement an 
operation HRz{a') with the wrong angle tan(a'/2) = (—1/3) * cot(a/2). To compensate this error, we measure the 
first qubit in the Z basis {\H), \V)} obtaining outcome ri and the second qubit in the basis B[{—1Y^ {a — a')] with 
outcome r2. When r2 = 0, which occurs with probability Ps(a — a'), we obtain the desired rotation HRz{a) up to a 
Pauli by-product operator. Thus, with such a block of two more qubits, we can boost the success probability of the 
rotation from Ps{a:) to Ps(a) ± [1 — Ps{o:)] -psia — a'). In the experiment, we have demonstrated such a building block 
to compensate the unwanted by-product operators induced by the wrong measurement outcome. In principle, one can 
generalize this technique to achieve a success probability arbitrary close to 1. Assume we have n blocks, the success 
probability will be lower bounded by Ps{a) + [1 — Ps(ck)][l — (1 — Pe)"] by repeating the above procedure. It can be 
verified that, for arbitrary e > 0, the extra computational overhead (the number of qubits) needed to achieve a success 
probability higher than 1 — e scales polynomially as 0(log i). The same analysis is valid also for the implementation 
of two-qubit entangling gates when we get the wrong measurement outcomes. Therefore, our experiment realizes a 
proof-of-principle demonstration of the trial-until-success strategy to compensate the randomness of the measurement 
outcomes in measurement-based quantum computation with non-cluster states. 



Additional experimental results of the two-qubit entangling gate 



In the experiment of the two-qubit entangling gate, when r2 = ^3 = 0, we use state tomography to characterize 
the output two-qubit state of the two-qubit gate. In Fig. |4j we list the experimental and theoretical density matrices 
when r4 = 1, Q = and r4 = 1, a = 7r/3. Also, we give the results when a takes the value of tt or 7r/2. 

When r2 7^ or ra ^ 0, we measurement qubit 4 in the Z basis {\H), \ V)} and decouple the two quantum wires. 
Rotation errors will be introduced in the operation, which can be compensated via the trial-until-success strategy as 
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demonstrated above. We have measured the fidchties of the experimental states and compare them with the expected 
ones, the results of which are listed in Table II when a is set to one of the angles {0, tt, 7r/2, 7r/3}. 



a a=^,r^ = \ b a = nil>,r^ = \ C a = 7xi2,r^ = ^ 





FIG. 4: The density matrices of the output states of the entangling two-qubit gate. In each sub-fig, the left-top: the real part of 
the experimental density matrices; the right-top: the real part of the expected density matrix; the left-bottom: the imaginary 
part of the experimental density matrices; the right-bottom: the imaginary part of the expected density matrix. The fidelity 
of each state is: a: 0.80 ± 0.03; b: 0.76 ± 0.04; c: 0.79 ± 0.03; d: 0.74 ± 0.04; e: 0.75 ± 0.02; f: 0.86 ± 0.02. 



DEUTSCH'S ALGORITHM 

In the main text, we consider only the case of ri = r2 = r^, = 7^4 = 0. Here we discuss in detail the cases when one 
or more outcome of four measurements is not 0. 

First, let us consider the implementation of the balanced function. We know that if we measure qubit 1 in 
the basis i?(7r) and qubits 2, 3, 4 in the basis i3(7r/2), balanced function can be implemented with the outcome 
ri = r2 = ra = 7-4 — 0. Below we discuss another two situations, (a) When r2 = ra = & (ri 7^ or r4 7^ 0), 
Pauli errors may occur on the output states, which will change the distribution of the final results. However, we 
can classically relabel the results to correct errors of the distribution. With regards to the classical relabeling, the 
detector signal corresponding to r2 = ra = 0, together with the measurements of the qubits 2,3,4 should be viewed 
as the entire function in the black box. In this case, the experimental results are shown in Fig. [S^ and the probability 
of successfully judging the types of the functions is 82% ± 3%. (b) When 7'2 7^ or rs 7^ 0, we can't implement the 
entangling gate with the given qubits and thus can't implement the balanced function. 

Next, we consider the case of implementing the constant function. When qubit 1 is measured in the basis Bij:) 
and qubits 2, 3, 4 in the basis S(0), we can simulate the constant function with the outcome r\ = r-i = ^ ^ 0. 
However, when the measurements are not zero, Pauli errors may occur on the output states and change the distribution. 
As with the implementation of the balanced function, we can also relabel the results. Considering again the case with 
the measurement results r2 = rj, = 0, in Fig. [SJd we show the experimental data when r\ and r4 take the value of 
either or 1. The success probability of recognizing the constant function is as large as 99% ± 1%. 
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TABLE III: The fidelities of the output states of the two-qubit entangling gate when we get the wrong measurement outcomes. 
We measure qubit 1 in the basis B\{a) to prepare the input logical state, and measure qubits 2, 3 and 4 in the bases 
S2(|-) = B-j,{^) = {s \H) ± ic |y) , c \H) — is \V)} and the Z basis {\H), \V)} respectively. The measurement outcomes are r2, 
rz and r^. The path degree of freedom of a photon is denoted as H' ^ 0, V' ^ 1 in the physical states. 




FIG. 5: a. The distribution of the final results when the function is balanced. The success probability of recognizing the 
constant function is 82% ± 3%. b. The distribution of the final results when the function is constant. The success probability 
of recognizing the constant function is 99% ±1%. 
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